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Summary 

■Will  upoit  piseciita  a  theoretical  analyais,for  the  behavior  of  long,  circu¬ 
lar,  cylindrical  shellfof  sandwich  construction  under  axial  compressive 
loads.  The  analysis  is  designed  to  evaluate  the  effects  of  the  relatively  low 
shearing  moduli  of  sandwich  cores  on  buckling  stresses.  Families  of  curves 
are  presented  for  use  in  designing  shells  of  sandwich  construction  having  Iso¬ 
tropic  facings  and  orthotropic  or  Isotropic  cores. 


The  results  of  the  theoretical  analysis  were  compared  with  those  obtained 
from  tests  on  a  series  of  curved  panels.  It  was  found  that  the  theory  applied 
reasonably  well  to  curved  plates  of  sises  sufficient  to  include  at  lenst  one 
ideal  bucUe.  Application  of  the  theory  thus  is  not  limited  to  long,  complete 
cylinders.  - 


^-This  progress  report  is  one  of  a  scries  prepared  and  distributed  by  the 
Forest  Products  Laboratory  under  V.  8.  Navy,  Bureau  of  Aeronautics 
Order  No.  NAsr  01237  and  01202,  and  17.  8.  Air  Force  No.  U8AF  18 
{&00)  -70.  Results  hers  reported  are  preliminary  and  may  be  revised 
as  additional  data  become  available.  Original  report  published  June  1932. 

•m 

■Maintained  at  Madieon,  Wie.,  in  cooperation  with  ibe  University  of  Wis¬ 
consin. 
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Introduction 


In  the  design  of  aircraft  am.  guided  missile*,  it  wai  found  necessary  to  da* 
vise  a  method  of  determining  the  etreaa  at  which  curved  aandwich  panels 
eubjected  to  axial  compression  become  elastically  unstable.  It  is  known 
that,  for  thin,  homogeneous  materials,  a  curved  form  greatly  increases  the 
critical  load  as  compared  to  a  flat  sheet  of  the  same  approximate  else.  A 
similar  increase  may  be  expected  for  curved  sandwich  panels.  Although 
this  report  applies  primarily  to  aandwich  construction  for  aircraft,  the  re¬ 
sults  are  general  and  apply  to  any  structures  of  the  type  considered. 

This  report  presents  a  theoretical  analysis  of  the  behavior  of  long,  circular, 
cylindrical  shells  of  sandwich  construction  under  axial  compressive  loads 
and  an  experimental  confirmation  of  this  analysis  by  teats  on  curved  panels 
of  sufficient  sise  to  Include  at  least  one  ideal  buckle.  Thus,  these  panels 
are  aatumed  to  simulate  the  action  of  complete  cylinder*. 

The  buckling  of  a  homogeneous,  isotropic,  thin-walled  cylinder  was  treated 
by  von  K*rman  and  Taien  (lb)—  and  by  Tsien  (IS.  14)  In  related  papers. 

These  authors  assumed,  in  addition  to  the  wave  form  of  the  classical  theory, 
inward  buckles  of  diamond  shape  to  represent  the  characteristic  buckles  that 
are  actually  observed.  They  used  an  energy  method  to  determine  the  criti¬ 
cal  compressive  stress.  This  method,  in  which  only  diamond- shaped  buckles 
are  user',  was  applied  by  March  (1)  to  cylinders  made  of  plywood,  an  ortho- 
tropic  material.  Particular  attention  was  paid  to  the  effect  of  initial  irregu¬ 
larities  that  contribute  to  the  obaerved  scatter  of  experimentally  determined 
critical  stresses  of  both  isotropic  and  orthotropic  cylinders. 

In  this  report,  the  effect  of  shear  deformation  in  the  core  of  a  sandwich 
cylinder  it  taken  into  account  by  employing  an  approximate  "tilting"  method.  - 
This  method  was  used  by  Williams,  Le,.gett,  and  Hopkins  in  their  analysis 
of  flat  aandwich  panels  (18)  and  by  Lej.,  vtt  and  Hopklna  in  their  analysis  of 
flat  sandwich  panels  and  cylinders  (4).  It  amounts  essentially  to  assuming 
that  the  transverse  components  of  shear  stress  are  constant  across  the  thick¬ 
ness  of  the  core.  The  form  of  buckles  assumed  by  Leggett  and  Hopkins  (4) 
in  the  cylinder  is  different  from  that  assumed  in  this  report. 

The  iors  and  facings  are  taken  to  be  orthotropic,  with  two  of  their  natural 
axes  parallel,  respectively,  to  the  axial  and  circumferential  directions  of 
the  cylinder.  The  facings,  which  may  be  tqual  or  unequal  in  thickness,  are 


—Underlined  numbers  in  parentheses  refer  to  Literature  Cited  at  end  of 
report. 
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iMumid  :o  be  thin,  but  th«lr  flexural  rigidities  are  not  neglected,  aa  theae 
may  be  of  Importance  in  certain  caaea.  All  atreea  componenta  in  *he  core 
are  neglected  except  the  transverse  ahear  componenta.  it  waa  pointed  out 
by  Reiaaher  (12)  that  the  atreaa  component  in  the  core  normal  to  the  facinga 
may  be  of  importance  in  the  analysis  of  sandwich  ahella.  Preliminary  calcu¬ 
lation*  Indicate  that  the  effect  of  thia  component  la  amall  in  the  problem  under 
consideration. 

As  was  done  in  work  described  by  Foreat  Products  Laboratory  Report  No. 
1322-A  (7),  initial  irregularities  are  assumed  to  be  present  and  to  grow  un¬ 
der  increasing  compressive  load  until  buckling  occurs.  For  a  discussion  of 
this  Important  matter,  the  reader  is  referred  to  that  report  and,  in  particu¬ 
lar,  to  the  observations  of  the  growth  of  artificially  produced  initial  irregu¬ 
larities.  Also  aa  described  in  report  No.  1322-A,  a  large  deflection  theory 
is  uaed  to  take  inso  account  the  nonlinear  support  associated  with  the  curva¬ 
ture  of  the  shell,  ae  discussed  by  von  Harman,  Dunn,  and  Tsien  (17).  The 
derivations  uf  the  differential  equation  for  a  stress  function  and  of  the  ex¬ 
pression  for  the  energy  of  deformation  are  extensions  of  the  analysis  used 
by  von  Kerman  and  Taien  (lb)  far  the  homogeneous,  isotropic  cylinder  to  the 
sandwich  cylinder  composed  of  orthotropic  materials.  Suitable  modification 
is  made  for  the  effect  of  shear  deformation  in  the  core  of  the  sandwich. 


Theoretical  Analysis 


Choice  of  Axes  Notation 


The  choice  of  axes  is  shown  in  figure  1,  the  coordinate  £  being  measured 
along  the  circumference.  The  notations  for  stress  and  strain  are  those  of 
Love's  treatise  {5).  The  components  of  displacement  in  the  axial,  circum- 
li.entisl,  and  radial  directions,  respectively,  are  u,  v,  and  w,  the  latter 
being  positive  inward.  Since  initial  irregularities  of  the  cylindrical  surface 
are  assumed,  the  symbol  wQ  U  used  to  denote  the  initial  distance,  measured 
radially,  of  a  point  of  the  middle  surface  from  a  true  cylindrical  surface  of 
radius  r,  and  the  symbol  w  to  denote  the  corresponding  distance  at  any  stage 
of  the  deformation.  The  thickness  of  ths  core  is  denoted  by  c  and  that  of  each 
facing  by  fj  and  f2,  respectively. 


Extensions!  Strains  and  Stresses 

Expressions  can  now  be  written  for  the  exteneional  strains  uniform  scross 
the  thickness  of  the  cylindrical  shell  and  for  the  corresponding  mean  mem¬ 
brane  stresses.  On  these  will  be  superposed  a  system  of  flexural  strains. 
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and  the  energy  of  dafertaiation  associated  with  each  ayatem  of  atralne  will 
be  found. 


The  extenslonal  atraine  era  expressed  by  the  equations i 
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In  each  facing,  the  corresponding  stress  components  are: 

E.. 


T  <* 


xx  +  °yx  eyy^ 


Ey 

Yy  a  IT  (eyy +  ff*y  e«) 


(2) 


X  =  u  e 
y  rxy  xy 


where  Ex  and  Ey  are  Young's  moduli,  pXy  is  the  modulus  of  rigidity  for 
shearing  strains  referred  to  the  x  and  £  directions,  a  and  a  are  Poisson's 

ratios,  and  \  ■  1  -  axy  ayx.  All  of  these  quantities  ere  elastic  properties  of 

the  faclnga.  Because  the  stress  components  X^,  Vy,  and  are  neglected 

in  the  core,  ‘he  mean  membrane  stress  components  for  the  cylinder  are* 

_  E. 

X,  .  (, 
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h  ■  c  +  fj  ♦  r2 


From  the  relation! 


Co  «  £  cr 
*  yx  y  xy 


(*) 

(6) 


that  hold*  for  orthotropic  material*  and  equations  (4),  it  follows  that: 

E  »  ■  E.  <r  (7) 

a  *yx  T>  xy 

By  uaing  this  relation,  it  is  found  from  equation  (3)  that: 
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(8) 


The  mean  membrane  streea  components  satisfy  the  equations  of  equilibrium: 


a3L  a* 

-r—  +  -r-^  *  0 

ax  ay 


ax  ay 


(9) 


They  can  consequently  be  expressed  in  terms  of  a  stress  function  as  follows: 

(10) 
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It  is  found  from  equations  (1)  by  eliminating  u  and  v^that: 
.2-  .2. 
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By  introducing  (10)  In  (0)  and  aubatitutlng  the  reaulta  In  (Id))  the  following 
differential  equation  for  F  la  obtained: 
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where: 


A  "4--  *  c  *—  --S2 
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03) 


It  la  readily  aetabliehed  that  the  following  expreaaion  repreaenta  the  energy 
of  extenaional  deformation  of  a  rectangular  portion  of  the  aheii  with  edgea  ot 

length  a  and  b; 
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dy  dx 


(M) 


Form  of  Bucklea  and  Initial  Irregularltlca 

The  street  components  XK,  Y  ,  end  X  In  (14)  are  derived  from  a  atreaa 
function  F,  aatiafying  the  differential  aquation  (12),  which  involvea  derive  - 
tivea  of  w0  and  w  representing  tha  inil<*l  ; r.i  d. formed  middle  aurfaca  of  tha 
■  hell,  TOY  w,  the  Inward  radial  deflection,  tha  following  form  will  be  choaan: 


—  ■  g  +  8  cos2  (fty  -  ox)  cos2  (0y  ♦  ex) 

T 

where 


(>s) 


0*) 


The  nodal  linea  of  the  trigonometric  portion  of  aquation  (15)  are  ahown  in 
figure  2.  The  diaplacament  w  is  positive  inward.  In  equation  (16),  a  and  b 
repreaent  the  length  and  widtli,  respectively,  ut  a  diamond.  The  initial  ir- 
reguiaritiea  will  be  aeeumed  to  have  the  form  (15).  Tbia  la  done  fos-  the  pur- 
poee  of  simplifying  the  calculationa,  Then  ia  choaan  in  the  following  form: 
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-jr  ■  s o  ♦  a0  co>2  (fir  -  <*)  co*2  (fir  *  «)  U7) 

An  initial  flat  spot  on  the  surface  of  the  cylinder  could  be  described  roughly 
by  equation  (17).  An  initial  irregularity  wae  introduced  here,  ae  it  wee  in 
report  No.  1322-A  (7),  to  obtain  a  qualitative  description  of  its  influence  in 
causing  an  iaolated  buckle  to  develop  In  ita  vicinity,  if  the  initial  depth  of 
an  irregularity  of  the  form  (17)  is  very  email,  the  dimension!  of  the  area 
that  it  occupies  are  not  very  important.  For  this  reason  in  order  to  simplify 
the  calculations,  the  dimension*  a  and  b  in  equation  (17)  are  taken  to  be  the 
same  as  those  in  equation  (IS),  From  the  qualitative  description  that  is  ob¬ 
tained  of  the  development  of  an  isolated  buckle,  conclusions  were  drawn  in 
report  No.  1322-A  that  led  to  the  derivation  of  the  final  formulas  from  the 
analysis  for  the  case  w0  «  0. 

Th*  details  of  substituting  (15)  and^(17]_in  {l2),__of  obtaining  the  stress  func  ■ 
tion  F  and  the  eircr*  component*  Xh,  V  ,  and  X,,,  of  substituting  these  stress 
components  in  equation  (14),  and  of  related  operations  are  identical  wish  the 
corresponding  operations  performed  In  report  No  1322-A  (7).  Reference 
is  therefore  made  to  equations  (21)  and  (31)  of  that  report.  The  following 
differences  in  notation  should  be  noted: 

Notation  of  Report  No,  1322-A  Notation  of  R-esent  Report 
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^yx  °xy 
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X,  Y,  Xy 

From  equation  (14),  the  energy  of  extcnaional  deformation  Wj  is  then  found 
to  he: 


W  ,.h*b 
W,  --g- 


jirhV  (S*  -  6„V 


J _  ,0  '  .¥  U  l 


128B01 


128Ao 


16  (Ao4  +  81B04  +  9Co*0*)  1 6  (8 1 A  <r4  +  B04  +  9Ctf*02) 

4  r iBi  (6  +  »oi)  I  2  2  8<r*y  1 

♦  — I - - 21 - U  4Bp  +  4Acj  *  Cj  ?  J 


64 (A a4  +  B04  ♦  C^) 
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This  equation  ie  equation  (31)  of  report  No.  1322-A  (7)  with  the  proper  veluee 
inserted  tor  the  former  abbreviations  M  and  8.  The  quantities  £a»d  Cj 

represent  the  mean  compressive  stress  and  the  mean  cireomferentiafTtress, 
respectively. 

If  n  is  the  number  of  buckles  in  a  circumference,  the  width  b  of  an  individual 
buckle  and  n  are  related  by  the  equation: 


b  .in 

(19) 

n 

Then: 

„  _  »  JL 

0*b  “  IF 

(20) 

It  will  be  convenient  to  denote  the  ratio  ~  of  the  dimensions  of  a  buckle  by: 


b  <t 

E  *  -  *  — 
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(U) 


in  the  expression  obtained  from  (id)  by  using  equations  (20)  and  (21),  let; 

2h  *  -  £  .  .  *  (22) 
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Equation  (18)  then  becomes: 
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(23) 

(24) 

(25) 

(26) 

(27) 

(28) 
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Flexural  Energy  of  the  Shell 


To  determine  flexural  energy  of  the  shell,  the  following  simplified  expres¬ 
sions  for  the  changes  in  curvature  and  unit  twist  are  used: 

•2  (w  -  *o)  »2  (w  -  W  )  s2  (w  -  w  ) 

w  w  v  ie  x\ 


A  discussion  of  the  approximations  Involved  will  be  found  in  a  paper  by 
Donnell  (1_),  These  expressions  were  used  by  von  Karman  and  'lslen  (16) 
and  by  March  (7).  The  expressions  (30)  are  exactly  those  used  in  calculat¬ 
ing  the  flexural  energy  of  a  flat  sandwich  plate.  The  approximate  flexural 
energy  of  such  e.  plate  was  found  by  March  (9)  and  by  Ericksen  and  March  (Z) 
by  using  the  "tilting"  method  of  Williams,  Leggett,  and  Hopkins  (4,  18). 


In  this  method  it  is  assumed  that  any  line  in  the  core  that  is  initially  straight 
and  normal  to  the  undeformed  plate  will  remain  straight  after  the  deforma¬ 
tion,  but  will  deviate  in  the  x  and  £  directions  from  the  normal  to  the  de¬ 
formed  plate  by  amounts  that  are  expressed  by  the  parameters  k  and  kL 
These  parameters  are  determined  by  an  energy  method.  These  "tilting" 
factors  k  and  k1  are  introduced  as  well  as  two  quantities  and  that  deter¬ 
mine  the  positions  of  the  surfaces  in  which,  respectively,  the  components  u 
and  v  of  the  displacement  in  the  core  vanish.  The  letters  k1  and  <£  replace 
h  and  £,  respectively,  of  report  No.  1583-B,  because  h  and  r  have  already 
been  uaed  In  the  present  report.  The  following  derivation  of  the  expression 
for  the  flexural  energy  follows  closely  that  used  for  the  flat  sandwich  panel 
in  Forest  Products  Laboratory  Report  No.  1S83-B  (2),  to  which  reference 
is  made  for  further  detail*.  For  the  sake  of  simplicity  in  writing,  the  initial 
Irregularity  wQ  will  be  for  the  present  taken  equal  to  sero.  It  will  then  be 

Introduced  in  the  final  atepa  by  replacing  w  by  w  -  wQ. 

The  components  of  displacement  in  the  core  (fig.  3)  are  taken  to  be: 

vc  k'  U  -  q’)~  (31) 

W£  a  w  (x,  y) 
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Thu*  l  *  q  denote*  the  surface  In  which  the  component*  of  displacement  In 
the  x  direction  vanish  and  h  It  the  parameter  describing  the  inclination  in 
the  x  direction  of  the  reepectlve  plane  section*  to  the  normal  to  the  deformed 
sur Ace.  Similarly  and  10  are  related  to  the  displacement*  in  the  ^  direc¬ 
tion.  These  four  quantities  are  to  be  determined  in  such  a  way  that  the  flexur¬ 
al  energy  associated  with  a  prescribed  deflection  win  minimum. 


To  arrive  at  expression*  for  the  components  of  displacement  In  the  facings, 
it  ie  noted  that  the  continuity  of  the  displacement  at  the  fecing-to-cora  bonds 
requires  that  the  componenta  {11),  evaluated  at  £  a  0  and  £  »  c,  shall  be 
those  at  the  inner  surfaces  of  the  facings  fj  and  f respectively.  Within 
each  facing,  the  components  of  dlaplacemlht  irOnumed  to  be  such  that 
a  straight  line  initially  normal  to  the  undeformed  surface  of  the  plate  will 
be  straight  and  normal  to  the  deformed  surface.  Accordingly,  the  compon¬ 
ents  of  displacement  in  the  facings  fj  and  respectively,  are; 

-  sw -ofg. 

gw 

vi 9  <k,ql  - 1> 

W,,  ■  w  (x,  y) 
and 

“2  ■  -[k  (c  ■  ’> 4 1  ■  'll? 

vi  ■-[•>■(*  -  q')  + 1  -  c]  g 

w2  a  w  (x,  y) 

The  componenta  of  strain  in  the  core  e  and  facings  f  j  and  will  be  denoted 
by  the  superacrlpts  c,  1,  and  Z,  respectively.  ~™ 


(32) 


(33) 


From  (31),  the  transverse  shear  strains  in  the  core  are: 


<i«l 


The  effect  of  the  remaining  strains  In  the  core  is  assumed  to  be  negligible. 

In  finding  the  strain  energy  of  the  facings  in  the  bending  of  the  plate  (or  shell), 
it  ie  convenient  to  consider  the  components  of  strain  in  the  facings  to  result 
from  the  superposition  of  two  states  of  strain-  The  first  of  these  consists  of 
the  membrane  strains  In  the  facing*  associated  with  flexure,  that  is  the 
strain  in  their  middle  surfaces.  From  (*2)  and  (JJ),  these  strain#  are  found 
to  be: 
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*xx  -  (k<«  ♦  y>;^r 

4>£ 

exy  "  ^  +  kV  +  V  |^7 


f2ld2w 


c  ■ 

XX 


k  (c '  q)  +  2  ax* 


<2)  f\  •  .  •  .  f/]  §2w 

%v  3- r  ,c-ql  +  T  “ 

L  J  91 

’»  ■  •  [mc  -  <0  *  k'  »C  -  ,')  *  ra] 


The  second  state  of  strain  in  the  facings  Is  that  associated  with  their  bending 
about  their  own  middle  surfaces.  This  statt,  in  either  facing,  has  the  com¬ 
ponents: 


'  _  .*  82w  »  _  i  92w  • 

e  -  ■  ci  c  —  -  C  ‘Vi  c 

**  8x2  yy  ay2  *y 


t'  tf?.  c'  =  -  21' 


where  £  is  measured  from  the  middle  surface  of  the  facing  under  considera¬ 
tion. 

The  strain  energy  in  the  core  or  facings  is  given  by  the  expression  (6,  8). 

U  m  TT  (f  f  E,  e2xx  +  E„  e2yy  4  2E„  see 
2XjfJj  *  )f  TI  X  yx  xx  yy 

4  Xp^  e2xy  4  >.pyt  e2yt  4  Xp^  dC  dy  dx  (38) 

where  for  the  material  under  consideration  (core  or  facing),  X  »  <  ?xyffyx: 
E„  and  are  Young’s  moduli,  p^y ,  py^,  and  p^x  are  moduli  m  rigidity; 

and  and  are  Poisson's  ratios,  Primetflelters  will  denote  the  elastic 

■y  y* 
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•*»  ? 


constant!  of  the  core  material  and  unprimed  letter!  will  denote  those  of  the 
facing  material.  The  Integration  indicated  in  formula  (38)  is  to  be  carried 
out  over  the  area  OABC  of  figure  2  and  the  thickness  of  the  core  or  facings. 

i 

The  energy  in  the  core  is  obtained  by  substituting  expressions  (34)  into  (38), 
the  remaining  strains  in  the  latter  formula  being  neglected  as  previously 
stated.  After  integrating  with  respect  to  l  over  the  thickness  of  the  core, 
the  expression  for  the  energy,  denoted  by  Uc,  is 

u=  '  f  J  j <>  •  M2  (£)'  ♦  ,\i  0  -  k’>2  (%)]  dv  d«  <W) 

The  etrsin  energy  in  the  facing*  *a&bela.ted  with  the  membrane  strains  is  the 
sum  of  the  energies  obtained  from  (35)  and  (3b).  With  the  substitution  of 
these  expressions  into  (38)  one  obtains,  after  integration  with  respect  to  t. 
the  foiiowing  expression  which  is  denoted  by  U w 


3  h  _  f 

"m-it/ /  E«V>(k’ 

n  n  *- 


*  Ey  fcV  *  'v  *  (k  •  V)  *  t)"}($9 ' 

*  2E«v{'i  (ki  *  7>(kV  -  y)  *  <2  (k  (c  -  O  *  y)(k  (C  -  q>  *  yjj^y^r 

+  k*‘xy|",l  *  kV  *  fj)2  *  <2  (k  (c  -  q)  t  k'  (e  -  <,')  ♦  dV  d*  (d0l 


The  strain  energy  in  the  facings  associated  with  the  flexural  strain,  Up,  is 
obtained  by  substituting  expressions  (37)  into  (38)  and  integrating  over  the 
volume  of  each  facing.  After  integrating  with  respect  to  (, 


vf-4^)// [■.(&/ -.gf/ 

*  fht’)  ]'■  ’■ 


ftic)  .it » •!»«?» 

Wi  * ^  ST  17 


dy  dx 
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Now  in  all  of  the  expressions  (39),  (40),  and  replace  w  by  w  -  w0.  The 
flexural  energy  of  the  region  OABC  of  the  shell  l§  the  eum  of  U^,  Up, 

and  U  •  “  ““ 

—  * 

For  equilibrium,  the  "tilting"  factors  k  and  kj_  and  the  ordinates  ^  and  of 
the  neutral  surfaces  are  to  be  chosen  so  that  the  total  energy  is  a  minimum. 
But  these  factors  appear  only  in  the  flexural  energy  Hence,  they  must 

be  chosen  to  satisfy  the  conditions  ~~~ 


By  proceeding  exactly  as  in  report  No.  1583-B  (£),  the  quadratic  form  (Al4) 
of  that  report  with  and  replacing  h  and  r,  respectively,  is  obtained  for 

Wjg.  The  coefficients  in  equation  (AI4)  are  dciined  by  equations  (Alb)  in 


termi  uf  the  quantities  A|,  which  are  defined  by; 
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fate-; 


On  factoring  out  the  common  factor,  it  la  found  that! 

2W2  -  r2  <6  -  60)2  ab£  Bj*  (kq)2  4  2B2'  (kq)  (k'q')  +  Bj'  (k'q')2 

+  2B4'  (kq)  k  +  28 g*  (kq)  k'  +  2B5'  (k’q  )  k  +  2Bfc'  (k'q')  k' 

+  B?'  k2  +  2B8'  kk'  4  B9'  k'2  4  2B1Q'  (kq)  +  2Bn'  (k'q) 

(47) 

where  the  quantities  Bj'  are  defined  in  terms  of  the  quantities  At*  by  equa¬ 
tions  (A15)  of  report  Mo.  1583-B,  each  At'  replacing  the  corresponding  A; 

in  those  equations.  The  quantity  in  brackets  in  equation  (47)  corresponds  to 
2U*  in  report  No.  1583-B.  (Note  that  equation  (A 22)  of  that  report  should 
read  P  =  2U*). 

It  is  easy  to  see  that  the  steps  of  imposing  the  conditions 


t  2Bj2  k  +  2B|j 


aw,  aw,  aw,  aw, 

a^kqT  a(k'q')  ek  ah' 


and  of  determining  kq,  k'q1,  k,  and  k^  and  substituting  their  values  in  the  ex¬ 
pression  (47)  for  2W2  are  identical  with  those  taken  in  report  No,  1583-B  (2) 

and  that  2W2  is  equal  to  the  right-hand  member  of  equation  (A25)  of  that  re¬ 
port  multiplied  by  r2  (6  -  6fl)2  ab.  It  is  concluded  from  equations  (A26),  (A27), 
and  (A28)  of  report  No.  1583-B  that: 

2  2  [  1  fAl'  *  2A2  +  A*  +  ^Al*A3*  * 

IX/  _  !  /  3  _2  /£  C  .  L  J  L  ”  i  ^*4  J 


W2  =  1/2  r4  (ft  -  b0)  ab 


•2wJL 

2  A . 


A.*  A  >  ♦2  A,  A,  -  A,'2) 

I  +  JL_  +  -I_  +  -  Vh - 

A4'  a5'  A  *A  ' 


+  I,  (A.*  +  2A2'  +  Aj') 


where  If,  and  £  are  defined  by: 
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(c  *  '»  *  ‘‘) 
’,*>!  \  1  1 


L.ii ul 

{  12 


c  ft  f2 


f!  +  f2 


and  Aj'  by: 

A,'  0  E.tt4  ♦  ^yffV) 

A,'  =  —  (E„ff.,w  +  \u„..)  o2ffZ 

£r  A  S  «  -  *  ~  * 

V  ■  i  (3  E/  * 


'  io2  f’tx  , 

A4  n - 8^’  A5 


8 


(49) 

(50) 

(51) 


(52) 


(53) 


(54) 

(55) 


Note  that 

A,'  +  2A  '  4  A  '  «  &-  K 
1  6  J  A  4 


where 


K4  ■  3E  r4  +  3E  4  2(E  «r  4  2Xp  )  t" 

4  x  q  '  x  yx  rxy' 

and  Introduce  the  following  abbreviations  in  the  expressions  for  Aj  ,  Ajj\ 


and  A^  : 

a  '  _  jSl  J  a'  fllj  A»  flj 

A1  X  dl-  A2  *  X  d2‘  A3  *  X  d3 


where: 


(56) 
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dJ  -  JE*  .4  +  X^y  *2 
h  "  <E*  *Y*  * 


d3  "  3Ey  +  ^Sty  * 

Alio, 

,  1  i  ,.2  i  2 

a'  .llliz .  ILlaL 

4  8  8 


(S?) 


(58) 


Substituting  these  expressions  for  A  j  equation  (48)  becomes 

4  ji  [k4  4  (dtd3  ■  d22)  er--^  + 

SS^di*  SS2^*  64fl4*2(d.d,  -  d,2) 

4  .  ■'■  .4.  .,  4  -  ■*-■-  4  -  *  _  A ... 

3*'yt 


‘LT  1  y  e  - 

A-S  n  2X  f“  ,6  '  6°'  abp 


*  w 


*  >  *  ?. 
%)*  Cx* 


?  f  t  J 


(59) 


The  fallowing  transformations  are  made  by  using  equations  (20)  and  (22): 

i£±- . 1  "2f _ .  _£at _ .JL* 

3kf't„  33“Vh  *.  “ 


where 


2EX* 

5»'n7^?h 


and 


where 

2E  ♦ 

s  - r“~ — 

y  3Xp rh 

rt 


(60) 


(61) 
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The  coefficient  of  the  expreaalon  in  bracket*  in  equation  (59)  ie  alao  trana- 
formed  by  using  equations  (20)  and  (22).  The  expression  for  W2  becomes i 

fivh5)  [m  ♦  (d,i,  -  d22)  ^  (|f  ♦  8^’ 


W 


(t  -  to) 
32kr2 


< 


1  ^i8x  ,  ^38y  ,  ^  (did3  "  d2Z).8x3y 


£  e‘ 
x 


E  2 
x 


♦  (Ig/h  )  Kj 


(62) 


Or 


W2  «  e4  ab  r\  (|  -  $„)2  \ 


(63) 


wncrc 


*4  =  32X 


1/»m-3\  i/jj 

j  i  \h  »  /  + 


1  + 


Tjd.S, 


1°X  r-i-y 

T?  ~ 


E. 


J  2,  JL  x  Q  tl 

‘ilE-(-*  VJ  ♦  (if/ h3j  k4; 


^x _ t; 


^,<13 

£  2  *2 


(64) 


Virtual  Work  of  the  Compressive  Load 


Exactly  aa  In  equation  (35)  of  report  No.  1322A  (7_),  the  virtual  work,  Wj, 
of  the  compressive  load,  calculated  for  the  region  OABC,  figure  2,  is 
found  in  the  notation  of  the  present  report  to  he: 


W, 


.  abh  [»p2  +  ^ar. 

=  abh  ffip2  + 


pci  +  IT  1-2  <6 1  ’  *o*>  p] 

pcj  +  e5  n  (t2  -  %0l)  p  7] 


(65) 


where 


e5 -htl 


(66) 


It  will  b  convenient  to  consider  the  mean  energy  per  unit  volume  of  the 
cyllndri  al  shell.  Hence: 


W  >  (W,  t  W?  -  W3)/abh 
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(67) 


In  accordance  with  aquations  (29),  (62),  and  (65)i 

w  ■  [•,-!*  (I2  -  t„2)2  -  «j0  (l2  ■  t02)  (t  -  t„)  ♦  «j  (t  -  t„)2 

+  ^n2  (t  - („>*]  ^  -  .Ji  tt2  -  e02)  p  7  m 


The  Buckling  Strega 

For  equilibrium,  the  derivative*  of  Vf  with  respect  to  the  varloue  parameters 
i>  ci>  6.-  and  a  vanish.  From  the  condition 

=  P,  it  follows  that  r,  *  ft  (69) 

8c  j  1 


Now.Cj  denote#  the  mean  circumferential  stress.  The  parameter  g  appears 

only  in  the  expression  for  Cj  as  given  by  equation  (30)  of  report  No.  13 22 -A 

(7).  The  fact  that  Cj  vanishes  implies  that  £,  which  describes  a  uniform 

radial  expansion  of  the  cylinder,  take*  on  such  a  value  that  the  mean  rircum- 
feren'JUl  stress  vanishes.  Further  consideration  of  the  parameter  £  is  not 
necessary. 


dW 

From  the  condition  ~~  *  0, 

*5 


it  follows  that: 


P 


*2  <H  +  €0> 

2ft 


,  "3  ,  •d’H  (ft  -  fto)  h 

*  4  J  'r 


(70) 


where  £,  as  previously  noted,  is  the  mean  compresalv*  stress. 
Let: 

*1  e2  *3  *4 

Vi  "{r,  Y,  -«r.  Y.  ■  «r  .  Y .  • -=r 


'1  Ea  '  t2  E,'  t3  E,  ’  t4  Et 

Then  (7  0)  can  be  written: 

„  T,  ,  A  t  ,  Y2  (3ft  ♦  ft0)  Yj  Ya*l 
p ■ «.  Uvp u  ♦  u  • — j|— 


({ •  t.)  ■< 

*s  7 


(71) 


(72) 


The  mean  compressive  strain  t_  is  expressed  by: 
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(7  3) 


A»  In  report  No.  1322»A  (?),  it  is  found  that: 

n  3  ^  2  *2  *  h 

*  ■  Bp  +  ^4  *  (t  »  t0  )  7 


The  further  equilibrium  condition* 


aw 

8*1 


•  0  end 


8W 

a* 


■  0 


ft  re  to  be  satisfied.  The  first  of  these  is  associated  with  the  number  of 
buckles  in  a  circumference  or  with  the  width  of  an  individual  buckle,  and  the 
second  with  the  -  alio  hit t  of  the  width  ^  buckle:  to  its  length.  It  i“  rid! 
analytically  feasible  to  use  these  conditions  in  connection  with  equation  (68). 


The  following  method  of  arriving  at  the  critical  value  of  £  is  based  upon  an 
extended  discussion  in  report  No.  1322- A  (7).  Briefly,  it  was  considered 
that  an  isolated  initial  irregularity  would  increase  in  size  and  depth  with  in¬ 
creasing  mean  compressive  stress  £.  It  was  iherefore  considered  that  the 
load-mean  compressive  strain  curve,  with  p  as  a  function  of  *^,  for  a  given 
small  initial  depth  of  irregularity  would  be  '.he  envelope  of  the  family  of 
curvet  for  i  function  of  £,  drawn  for  i  series  of  values  of  tj  by  combin¬ 
ing  equations  (71)  *nd  (72),  On  taking  into  consideration  the  possibility  of 
jumps  from  one  energy  level  to  another,  i!  was  concluded  that  the  critical 
values  of  £  tvocld  scatter  considerably,  as  they  actually  do  in  test,  depend¬ 
ing  upon  the  depth  of  the  initial  irregularity  and  the  characteristics  of  the 
loading  process.  It  was  noted  that  the  value  of  £  at  the  relative  minimum 
point  on  the  envelope  of  the  curves  for  £  as  a  function  of  t ,  drawn  for  £0  n  0, 
was  intermediate  among  the  possible  critical  values  of  £.  This  mlnirnuhn 
was  accordingly  chosen  as  the  "theoretical”  critical  stress,  because  it  could 
be  conveniently  determined  by  finding  s  relative  minimum  of  £  as  a  function 
of  £  and  rj.  It  is  necessary  to  employ  numerical  methods  to  determine  the 
relative  minimum  value  of  £. 

In  report  No.  1322- A  (7),  the  aspect  ratio  t  of  the  buckles  was  assumed  on 
the  basis  of  experimental  observations  before  the  minimisation  of  £  was 
undertaken,  Here,  because  of  the  influence  of  shear  deformation  in  the  core, 
a  suitable  value  to  assign  to  s  can  not  be  estimated. 

Equation  (71),  with  a  0,  can  be  written  !n  the  form; 

P  »  XE,  7  (~ *t) 
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where 


k-^iK  -  v  -r  *  ^4 "] 


The  mean  compreaaive  atreaa,  pf,  in  the  facing*  ia  related  to  the  mean  com- 
preaalve  atreaa  in  the  ahell  by  tKl  equation 

Pf  (fl  +  *2) 

p  a  -I — 1- -  (76) 

On  recalling  the  definition  of  Ea,  it  ia  seen  that  equation  (74)  can  be  written 
PfaKEK£  (77) 

For  a  relative  minimum  of  p^,  the  condition 
■?—  ■  0  muit  he  latiafied. 

ee 

From  this  condition,  it  follows  that; 


The  substitution  of  thia  value  of  £  In  (75)  ylelda 
64  ,V3  9>22 

K'£i<-ir5F^r  ♦'«'>>  <7,) 

In  equation  (79),  K  ia  a  function  of  n  and  a,  which  occur  in  the  definitions  of 
3,  Yj,  y^,  y3>  and  By  using  the  definitions  of  the  quantities  E^,  E^, 

and  pm~tEat  appear  through  the  symbols  A,  B,  and  C  in  the  equations  f2&)i 
(27),  and  (28),  the  following  expressions  are  obtained  for  Yj<  V^>  and  yj 
(see  equations  (71),  (26),  (27),  and  (28):  ~  " 


_ +  _ J_ _ , 

v]  "4096  4096  Ex 


4096  4096  E„  .  E-  o.2E  , 

*  512  (s4=£  +  81  18*  z2) 

y  Pxy  y 

+ _  »4 _ + _  n«4 _ 

4  Ex  ,  Ex  ,  “Ejt  72eI  7 

512  (Bit’  £-+  1  +  9**s-  -  lBtr  a*)  2048  (a4  =-  4  1  +  ~*-2r  tT) 

y  y  1  y  *xy  y 
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(81) 


1  +  a 


17  a4 


1  4096  512  (a2 +  9)2  512  (9a2  +  l)2  2048  (1  +  a2)2 

.4 


1 


2  512  32  (1  +  *2)2 


*3  *  256  32  (1  +  *2)i  .  , ,  .  .2.  e 


4  in  ».2 


32  \  h*  {fj  +  C2) 


(3a4  4  3  4  2s2)  E2 


ndjSx  T*3Sy  t^Mv  (dld3  -  d22) 

,  + - —  ♦  - *■  +  - **T  '7 - 

E  fc*- 


Ea‘ 


(87) 

(88) 

(89) 

(90) 


where 

T  =  3z4  t  3  4  2z2 


a 


2  $  F. 

3X^*rh 


2  6  E 

sy *  3>^Vh 


(91) 

(92) 


Alter  icme  manipulation  involving  aubatitution  of  expreaaiona  (or  yj,  y^, 

\y  and  y^,  formula  (79)  for  K  for  aandwtch  construction  with  iaotropic 
facinga  and  orthotropic  core  can  be  written  aa; 


M1  21 
•»  i  _ 

M2,  +  M3r,2Sx  f  If  „  1 

(93) 

^  2  _  ^  > 

11  3X  "  (fj  *  f2) 

.1  +  m4t)Sx  +  M5n2sx2  1  J 

where 

64  ,  9Y22 

M1  8  3x2  (v3  '  32Yj  ' 

(94) 

M2  ■  -J 

(95) 

E2  a*  ^ 

(96) 
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dl  1 

M4  >  (-J+  d3«)  2 

(97) 

<Vj  -  »2*>  ° 

M5  "  e2  ,2 

(98) 

Sy  H  tx 

0  ■  n1  or  0  ■  -~- 

S*  H  yt 

(99) 

dj  ■  3Es4  ♦  j  (i  -  <r)  £*2 

(100) 

d2  *  E«rs2  +  ~  (1  -  a)  Es2 

(101) 

dj»3E+~ (1-cr)  Es2 

te 

(102) 

a 

If  v  Is  taken  to  be  — ,  then: 

*  3Es2  (s2  +  ^) 

(103) 

ii  ■  s  E,i 

(104) 

'S  ' !E  <T*  11 

(105) 

If  tlio  V 1  •  Y^I  and  £  are  expressed  in  terms  of  s  and  0  (Eq.  87,  88,  89, 

and  90)  then  the  following  expressions  can  be  used  in  formula  89: 

3  1  2 

<106) 


3  _L  +  « 

U*-.  4(1  ♦»*)* 


M,  .  '  - - jft-j liL±»XI _ 

1  12s*  3(1  +  s2)2  1  +  tA  4  s4  ^  a4  4 


17i. 


M2-3.^2+~ 


128  k6(s24  9)2  16(9»2  +  l)2  64(1  +  *2)2 

(107) 


M3  -  ~  (9s4  a  70s2  +  9)  (p-  +  0) 


(108) 
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M 


M4  .  |  [%Ml  +  1  ♦  (»Z  ♦  ®)  *] 

M5  ■  5  [9t4  +  70* 2  +  9]  e 


(109) 


(no) 


For  construction*  for  which  the  shear  deformation  in  the  core  i«  negligible, 

f.  1.  1.  wh.n  U  v.„  l.r„  ,nd  i  I.  «»«..  »,  -n.y  b.  »k.«  .,«bl  ..  ..... 
■»* 

Then  expression  (93)  can  be  minimUed  with  reepect  to  j>  re*ultlng  in: 


K0«  2 


(1+1^ 

3  k  (fj  +  f2)  h2 


Thue,  K0  i*  a  function  of  end  M*  end  the  .tiffne.e  of  the  eendwich.  It 
wae  found  by  computation  that  a  relative  minimum  of  MjM2  =  0,24  occurs  at 
*  *  0.95.  The  minimum  buckling  -tree*  ia  then  proportional  to: 


°  5Q, 

where 


/Mfi  ♦  h) h  (113) 

By  letting  N  =»  jr-,  the  following  expression  can  be  written  from  equation  (89) 
for  constructions  having  any  value  of  Sxt 


SMjQj  5 

N  »  — TT +Trs' "71 


.  u  -a  j.  li  -«a  ^ 


1<1+C*2'  1  +  M4t*k  *  MjtlV 


where 

k  (11! 

°2  "  I 

II  wa*  found  in  Foreat  Products  Laboratory  P-eport  No.  1505  (10)  that  the 
values  of  I  t  lj  and  lj  can  be  expressed  as  follows: 
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(116) 

(11?) 

where* 

(118) 

After  substituting  these  expressions  in  the  formulas  for  Qj  end  Q£  end 
simplifying  Qj  and  become: 


(1  -  c3/h5)  (I  -  c/h)  - 


Ucd 


(U9) 


i  (1  -  c/h)4  +  3  (l  -  c/hfp 

- - -j- 

(1  -  c3/h3)  (1  -  c/h)  -  ^  (1  .  C/h)4  -  3(1  -  c/h)*  "  - 


(120) 


In  equation  (114)  lw4^>  M3.  M^,  end  Mj  depend  upon  s  and  U,  end 

end  Q^.  depend  upon  c/h  and  d/h.  Formula  (114)  can  Ihen  be  written  with 

appropriate  values  of  8,  c/h,  and  d/h  and  then  a  relative  minimum  value  14 
found  by  choosing  a  series  of  values  of  a  and  ij.  The  facing  stress  at  which 
buckling  will  occur  is  then  given  by: 


(121) 


Buckling  Stress  of  Sandwich  Constructions 
with  Isotropic  Facings  of  Equal  Thickness 
and  Orthotropic  or  Isotropic  Core 

Factors  in  formula  (114)  can  be  simplified  for  sandwich  constructions  hav¬ 
ing  facings  of  equal  thickne....  *..«n  d  ■  0  and  after  simplification 


Q, 


+  i  lh  +  1 


(122) 
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*y  >vj,  v^-AV  V  fc^^r*K37?yMrvir- 


I  /i  -  c/h\ 2 

1  \TT77hy 


and  finally  formula  (114)  become# 


5Qt  )m1  (1  +  c/h)2  M2^  *  MlAt 

N  "  4  n  24K  1  +  M  qS  ♦  M,q2 

4  x  5 


♦  (‘  Tir1'  “il  <1M> 


The  buckling  load,  which  it  proportioned  to  N  is  obtained  by  finding  the  lowest 
relative  minimum  of  expression  (124)  with  respect  to  t]  and  s.  Expression 
(124)  can  be  minimised  by  taking  a  derivative  with  respect  to  jrj  and  setting 
the  derivative  equal  to  sero.  This  leads  to  a  sixth  power  equation  in  q. 
Minimum  root!  of  jj  with  respect  to  a_  and  for  various  values  of  S*,  c/h,  and 
S  were  determined  by  means  ef  a  digital  computer.  Minimum  vain* a  of  N 
at  varloua  valuea  of  S*  and  for  c/h  equal  to  0.9,  0.8,  0.7,  and  0  equal  to 
0.4,  1.0,  and  2.5  are  given  in  Table  1  and  shown  as  functions  of  Sx  in  figures 

4,  5,  and  6.  Also  included  in  the  table  and  figures  are  values  of  N  for 
c/h  =  1.  These  values  represent  sandwich  constructions  for  which  the  stiff- 
ness  of  the  Individual  facings  are  aasumed  to  be  zero.  Although  no  actual 
constructions  can  be  made  of  thia  type,  the  values  can  be  considered  as 
representing  the  limit  for  constructions  having  extremely  thin  facings.  These 
values  of  N^were  obtained  as  follows.  Substitution  of  c/h  c  1  in  equation  (124) 
for  N  leads  to 

r. .  .  . .  ...  2.  v  ^ 


5  r— 

n.jVT 


1  /  M^q  +  Mjh 
V  ♦  m,tSx  t : 


M5q2Sx2> 


which  has  one  relative  minimum  value  for  q  ■  «x  This  minimum  value  is 
given  by 

5Mi 


Substituting  in  this  squation  the  valuea  of  Mj  and  Mg  given  by  equations  (108) 
and  (110)  yields  “ 

5  (~Z  +  •) 

N- — - -  (127) 

1  i/T  8X  9 
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which  is  minimum  for  ■  ■  «a  This  minimum  value  (for  »  ■  i/4;  hence 
k  ■  15/16)  i*  given  by 


3{rTax 


(128) 


Subetitution  of  thie  value  of  N.  in  equation  (1 18)  and  ueing  the  value  of  8X  for 
fj  •  fj  and  c/h  »  1  leads  to  the  following  limiting  expression  for  p^i 


4£h  5  n  fi*'cx  _  h 

p«  *  '  Ip-' u 


(129) 


For  values  of  S K  ranging  from  0  to  about  0<6  it  was  found  that  equation  (128) 
did  nnt  gi«»  InSSii  minimum  values.  In  this  ranee  of  the  minimum  values 

were  obtained  from  equation  (124)  by  use  of  a  digital  computer. 

The  value  g*  H  given  by  equation  028)  and  the  value  of  the  stress  given  by 
equation  (129y~are  independent  of  the  radius  of  the  cylinder  and  are  the  usual 
critical  values  associated  with  shear  Instability  of  the  core  (15). 

The  value  of  8  of  0.4  and  ita  reciprocal  2.5  were  used  in  the  calculations  be¬ 
cause  they  apply  to  honeycomb  cores  oriented  with  the  weak  direction  and 
the  strong  direction  parallel  to  the  length  of  the  cylinder.  It  has  been  noted 
from  figures  4  and  6  that  in  the  range  of  small  values  of  SK  where  N  is  inde¬ 
pendent  of  c/h,  the  orientation  of  the  core  makes  little  difference  m  the 
value  of  N. 


Application  of  Theoretical  Results 

The  compressive  facing  stress  st  which  buckling  of  cylinders  of  orthrotropic 
ssndwich  conetruction  occurs  is  given  by  equation  (77). 

P("KE»r 


where  Eu  ie  modulus  of  elasticity  of  facings  in  axial  direction,  h  is  sand¬ 
wich  thickness,  r  is  mean  radius  of  curvature,  and  Kit  given  by  formula 
(79)  »• 


K 


-*±Fii 

n 


*v22 

^ *  Vi 
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whsre  values  of  Yg  functions  of  _a  accoiding  to  aquations  (80),  (81), 

and  (82)  and  ia  a  function  of  tj  and  a  according  to  aquation  (86)  and  K  it 
taken  aa  tha  laaat  relative  minimum  with  respect  to  ij  and  *. 


For  sandwich  constructions  having  isotropic  facings  of  unequal  thickness 
and  orthotropic  cors  K  is  given  by 


X 


4N 

5Q1 


where  N  is  given  by  aquation  (114)  as 


iUxQv  % 

4rj'  +  60^(1  +  02) 


1  +  M4qSK  +  MjtfS*2 


+  Q2M2v> 


where  Qj  and  are  given  by  equations  (119)  and  (120)  and  Mj,  My 
M4,  and  are  functions  of  £  according  to  equations  (106),  "fib?}"  (108)', 
(T09),  andTllO)  and  N  is  taken  aa  the  least  relative  minimum  with  respect 
to  t|  and  £. 


For  sandwich  constructions  having  isotropic  facings  (Poisson's  ratio  1/4)  of 
equal  thickness  and  orthotropic  core  such  that  6  *  O.iior  8*>2,5i  or  isotropic 
core  (9  >  1.0)  equation  (1 14)  for  N  has  bean  solved  for  t/h  =  1,0,  0.9,  0.8, 
and  0.7.  Values  of  N  for  various  Sx  values  are  given  in  Table  1  and  in 
graphs  in  figures  4,  3,  and  6.  Then  the  critical  facing  stress  is  given  by 


4N  h 
rt  “  507  E  r 


where 


Q 


1 


3  /  5 

2yc"27h2  +  c/h  + 1 


and  N  is  given  In  terms  of  Sx  where 


—These  ratios  for  8  were  chosen  as  representstlv*  of  honeycomb  cores  such 
as  were  evaluated  in  Forest  Products  Laboratory  Report  No.  1849. 
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UcfC 


4Srh*' 


W 


and  £  Is  core  thickness,  {  is  teeing  thickness,  £  is  modulus  of  elestlcity  of 
f#;in^a,  h  is  sandwich  thickness,  (c  4  2f),  r  is  mean  radius  of  curvature, 
and  |i  ,  is  modulus  of  rigidity  of  core  associated  with  shear  strains  in  the 
_ k* 

axial-radial  plane. 


The  graphs  can  be  used  with  little  error  for  determining  N  for  constructions 
having  facings  of  unsqual  thickness,  provided  Sx  is  calculated  using  formula 

(60)  and  Qj  is  calculated  using  equation  (119). 

The  analysis  may  be  extended  to  apply  at  stresses  greater  than  the  propor¬ 
tional  limit  stress  of  the  facings  by  use  of  an  appropriate  tangent  or  reduced 
moduJu*  of  elasticity  for  the  facings.  This  entails  a  "trial •=  and- error"  solu¬ 
tion  involving  use  of  the  tangent  or  reduced  modulus  in  the  quantity  S  and 
elsewhere  until  the  resultant  facing  stress  is  compatible  with  the  stress- 
modulti*  Curve. 


Results  of  the  theoretical  analysis  fall  approximately  into  three  rones,  de¬ 
pending  upon  whether  there  is  no  shear  deformation  in  the  core  (Sx  *  0), 
some  shear  deformation  in  the  core  (small  values  of  S*),  or  considerable 
shear  deformation  in  the  core  (large  values  of  Sx). 

For  no  shear  deformation,  the  buckling  stress  is  determined  essentially  by 
means  of  the  isotropic  or  orthotropic  theory  (depending  upon  facing  proper¬ 
ties)  with  the  stiffness  determined  by  considering  the  spaced  facings  of  the 
sandwich. 


For  large  shear  deformations,  the  critical  stress  is  associated  with  insta¬ 
bility  of  the  core  in  shear.  This  hss  been  obssrved  for  sandwich  construc¬ 
tions  in  gensral  (IS),  and  it  has  besn  found  that  the  mean  critical  atress 
thus  determined  is  the  same,  regardless  of  the  original  assumption  of  the 
buckled  shape.  The  amalleat  value  of  S^at  which  the  critical  streas  is 

determined  by  shear  instability  of  tht.  core,  however,  is  greatly  affected  by 
the  assumed  form  of  the  buckled  shape.  The  inclusion  of  the  stiffnesses  of 
the  facings  If  gives  rise  to  the  family  of  curves  for  different  values  of  c/h, 

as  shown  in  figures  4,  5,  and  6,  Instead  of  a  single  curve.  If  the  atlffneeses 
of  the  individual  facings  had  besn  neglected,  one  curve  only,  that  for  c/h  ■  1, 
would  have  resulted.  The  percentage  Increase  in  buckling  stress  due  to  the 
stiffnesses  of  the  individual  facings  increases  as  the  shear  deformation  in¬ 
creases.  For  small  shea,  deformations,  the  increase  is  negligible. 
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From  the  rtUMtlng  involved  in  the  theory  leiMng  to  equation  (77),  consider¬ 
able  ecetter  tn  the  experimental  values  of  tints  ritloal  stress  is  to  he  expected 
because  of  the  effect  of  initial  irregularities.  fiSlmilar  scatter  Is  exhibited 
by  homogeneous  cylindrical  shells,  for  the  sweat  reason. 

The  possibility  of  failure  by  wrinkling  of  the  ti'^cinge  at  a  stress  lower  than 
that  predicted  by  equation  (77)  should  be  considered. 

The  analysis  in  this  report  involves  a  numbsiooi  approximations  and  assump¬ 
tions.  Such  procedures  are  necessary  until  tenors  rigorous  treatment  of  the 
problem  is  developed.  A  completely  rtgoroui  treatment  of  the  buckling  of  a 
homogeneous  cylindrical  shell  is  still  lacking,  In  spite  of  the  noteworthy 
contributions  of  von  Karmen  and  Tsien. 


Tests  of  Carved  PsMfcls 


The  large  site  of  complete  circular,  cyundric  ai  sheila  having  realistic  !sc- 
ing  and  core  thicknesses  inti  curvatures  could  not  be  adapted  to  the  available 
testing  apparatus,  Therefore,  axial  compreiswive  tests  were  conducted  on 
rectangular  panels  curved  to  various  radii,  F"he  dimension!  of  these  panels 
were  chosen  so  that  their  widths  and  lengths  *~*ere  large  enough  to  include  at 

least  one  buckle  of  a  else  predicted  by  theory  and  a»~-- )  as  shown 

in  table  2.  <  ■  n 

it  was  then  assumed  that  the  curved  panel  wou  ld  behave  approximately  as  a 
complete  cylinder.  The  type  of  edge  support  (^described  later)  was  such  as 
to  produce  no  clamping. 


Test  Specimen# 

The  teat  specimens  were  essentially  of  isotio  pic  construction  having  facings 
of  clad  24ST  aluminum  alloy  on  cores  of  eitht  t  balsa  wood,  oriented  so  that 
the  grain  direction  was  normal  to  tha  facing!,  or  of  corkboard  of  three  dif¬ 
ferent  densities.  Corkboard  cores  were  ch«l  an,  because  their  low  moduli 
of  rigidity  afforded  means  of  exploring  shelli  in  which  sis  cable  reductions 
of  buckling  stresses,  caused  by  large  core  i)K.esr  deformations,  could  easily 
be  obtained.  These  corkboard  cores  had  shumring  moduli  of  1,500,  950,  and 
520  pounds  per  square  inch,  as  compared  to K.  5,000  pounds  per  square  inch 
for  the  end-grain,  balsa-wood  core. 
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Dimensions  of  tha  specimens  arc  (Ivan  tn  table  2.  The  panel  aiaea  ranged 
from  approximately  70  inchca  square  to  panela  12  inchea  wide  and  30  inchea 
long,  Mean  radii  of  curvature  ranged  from  approximately  90  inchea  to  10 
inchea.  The  aandwich  c<matructlona  had  facings  of  0.012  inch,  0,020  Inch, 
or  0.032  Inch  thickness  on  cores  of  approximately  1/8  inch,  1/4  inch,  or  1/ 
inch  thickneaa.  All  constructions  tested  had  facings  of  equal  thickness. 

The  specimens  were  manufactured  by  the  bag-molding  process.  Detailed 
description  of  techniques  and  bonding  adhesives  used  in  this  process  are 
given  in  Forest  Products  Laboratory  Report  Ns.  1574(3).  The  curvature 
wae  attained  at  the  time  of  molding  by  uelng  a  steel  mold  curved  to  the  de- 
atred  radius.  A  strip  of  aluminum  1  inch  wide  and  0.032  inch  thick  was 
bonded  to  the  facinge  at  each  end  of  the  specimen.  This  was  done  to  facili¬ 
tate  machining  of  the  epecimen  ends  and  aleo  to  prevent  local  end  failure 
during  the  teat.  The  ends  of  the  specimens  were  machined  square  and  true 
in  a  milling  machine. 


Teeting 

The  vertical  edges  of  the  specimens  were  held  straight  by  loose-litting  wood 
guides.  These  guides  were  approximately  2  inches  by  2  inches  in  cross  sec¬ 
tion  and  of  lengths  1/4  inch  shorter  than  the  test  specimen.  They  were 
grooved  in  the  lengthwise  direction  with  grooves  approximately  1/4  inch  deep 
and  wide  enough  to  allow  the  guides  to  be  slipped  onto  the  edges  of  the  teat 
specimen.  No  attempt  was  made  to  clamp  the  vertical  edges  by  fitting  the 
guides  tightly. 

The  lower  ends  of  specimens  not  wider  than  30  inches  were  placed  on  a 
heavy  flat  plate,  which  was  supported  by  a  spherical  bearing  placed  on  the 
lower  head  of  a  hydraulic  teeting  machine.  The  heads  of  the  testing  machine 
were  then  brought  together  until  the  specimen  just  touched  the  upper  platen 
with  no  load  indicated.  Adjustments  were  nr  ao.j  on  the  spherical  base  until 
no  light  could  hr  seen  between  the  ends  of  the  specimen  and  the  loading  heads, 
Screw  jacks  were  then  placed  under  the  lower  loading  plate  to  prevent  tilting 
of  the  plate  while  the  load  was  being  applied  to  the  specimen.  A  single  thick¬ 
ness  of  blotting  paper  was  Inserted  at  the  ends  of  the  specimen  to  help  pre¬ 
vent  Iocel  end  failures.  The  load  was  then  applied  slowly  until  failure  oc¬ 
curred. 

Specimens  wider  than  30  inchea  were  tested  between  the  heads  of  a  four- 
•  crew,  mechanically  operated,  testing  machine.  No  spherical  bearing  was 
used.  The  specimens  were  cut  as  true  ae  possible.  If  light  could  be  seen 
between  the  ends  of  the  specimen  snd  the  heads  of  the  testing  machine,  shims 
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of  paper  or  brtu  w.eru  Inserted  until  the  gap  waa  closed.  That*  wide  sptei- 
mana  were  also  vary  loaf)  thsr  store,  a  mall  irregularities  la  the  and  bearing 
wara  abserbsd  early  In  the  test  without  causing  largo  variation#  from  uni¬ 
formity  In  the  ctraaaea  In  the  facing  a. 


Result#  of  Teat# 


The  facing  stresses  at  the  failing  loads  of  th*  curved  panala  are  given  in 
table  2.  For  later  comparison  with  theoretical  value#,  the  parameter  N 
was  calculated  for  each  teat  epeclm«n  by  using  the  formula  *" 


N  * 


5QlPfg 

4&ti 


where 

E  b  10.000.000  pounds  per  square  Inch  (modulus  of  elasticity  of  facings) 

The  visible  failures  of  the  specimens  were  of  a  type  caused  by  buckling. 
Large,  thin  specimens  actually  showed  large  buckles,  which  disappeared 
after  release  of  load.  The  appearance  of  these  buckles  always  caused  a 
sudden  drop  In  the  load.  Small,  thick  specimens  showed  buckling,  followed 
immediately  by  a  crimping  appearance  at  the  edges  of  the  buckle.  This 
crimping  was  undoubtedly  due  to  shear  failure  of  the  core  caused  by  high 
stresses  induced  in  the  sandwich  by  the  buckle.  Many  of  the  thick  speci¬ 
mens  exhibited  no  visible  signs  of  buckling  but  showed  similar  crimping. 
The  rapidity  of  failure  occurring  immediately  upon  buckling  undoubtedly 
prevented  visual  observation  of  the  buckle  itself.  Similar  behavior  was  ob¬ 
served  for  cylindrical  shells  of  plywood  (II). 


Comparison  of  Theoretical  and  Experimental  Re  suite 

The  theoretical  and  experimental  values  of  N  ere  given  in  table  2.  A  com¬ 
parison  between  them  may  be  obtained  by  referring  to  figure  7  which  shows 
the  experimental  values  plotted  against  the  theoretical  values,  The  scatter 
of  points  about  a  line  representing  equality  between  experimental  and  theoretl 
cal  values  ehows  that  theory  and  experiment  agree  within  approximately  +  30 
percent, 

In  view  of  the  inevitable  scatter  of  experimentally  determined  buckling 
■  tresses  that  is  associated  with  initial  irregularities  of  shape  and  variations 
of  material  properties,  it  is  concluded  that  the  agreement  between  results  of 
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teats  and  of  theory  1(3  satisfactory.  The  scatter  i«  no*  as  great  for  thetts  of 
iwdwibh  cMitmtlon  ae  that  okssrved  for  thin,  homogeneous  shells  $lg.  44, 
report  No.  IU2>A  ft))  or  for  plywood  shells  (fig.  3,  Forest  Products  Labors 
tory  Report  No.  132? (ID).  This  redaction  In  scatter  may  ba  attributed  W  a 
greater  total  thickness of  shall.  That,  irregularities  that  dspart  from  the 
true  cylindrical  surface  of  the  erdar  of  the  thickness  of  the  shall  art  less 
likely  to  occur  in  eandwieh  shells  than  in  thin,  homogsntous  shells. 


Conclusions 


The  buckling  stress  of  long,  thin-walled,  circular  cylinders  of  sandwich  con¬ 
struction  in  axial  compression  can  be  found  with  satisfactory  accuracy  by  the 
formulas  and  curves  of  tht  approximate  theoretical  analysis  of  this  report. 

Curved  panels  of  ei»€ii  large  enough  to  inrinde  at  least  one  ideal  buckle 
..  itr  Zwr .  , 

is  38  “jp  »sd  a  buckle  at  stresses  approximately  equal  to  those  of  a 

long,  complete  cylinder. 
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Notation 


A 


S 
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length  of  buckle. 

1/Eb. 

defined  by  equations  (42)  to  (46). 
defined  by  equation*  (52)  to  (55). 
width  of  buckle. 

1  /eb. 

thicknesa  of  the  core, 
mean  circumferential  §tr*«*. 

1  2<rxy 
ilrn  ®a 

fl  *  f2 
2 

defined  by  equation  (57). 
components  of  strain. 

defined  by  equation*  (26),  (27),  and  (28). 

defined  by  equation  (64). 

defined  by  equation  (66). 

Young’*  moduli  of  the  facing* 

Ex  <<1  ♦  ll ) 
h 

Ey  (f,  +  f2> 
h 

thickne»*e*  of  the  facing*. 

quantity  proportional  to  mean  radial  expansion. 
-56- 


h 

1 


k,  k* 
K 


U„  U, 


tt. 

dtfintd  by  aquation  (49). 
defined  by  aquation  (40). 
parameter  a  introduced  in  equation*  (31). 
lee  equation*  (74),  (75),  and  (79). 
defined  by  equation*  (23),  (24),  and  (25). 
defined  by  equation  (84). 
r/b  . 

mean  compreacive  stress, 
compreeaive  stress  in  the  facing*, 
introduced  in  equation*  (31). 

radiu*  of  middle  surface  of  the  cylindrical  shell. 

defined  by  equation*  (92). 

defined  by  equation  (83). 

axial  component  of  displacement. 

circumferential  component  of  displacement, 

•  train  energy  of  the  core  in  the  bending  of  the 
aandwich  (hell. 

■train  energy  of  the  facing*  in  the  bending  of  the 
aandwich  ahell. 

redid  component  of  displacement, 
extensional  strain  energy, 
flexural  strain  energy. 


# 


i 

l 


i 


i 
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a  -  ^ 


Xyl  Xy,  «tC. 


virtual  work  of  tho  eompraaaiv*  load. 
(Wj  +  W2  -  Wj)/abh. 
component*  of  etreee. 


V  V  y3.  y4 


K  tx’  H  yt 


defined  by  equation*  (80),  (8 1  )*  (82),  and  (86). 

»  parameter  that  ia  proportional  to  depth  of  a 

buckle. 

initial  value  of  6.  > 

mean  compreeeive  etrain. 
coordinate  ihown  in  fig.  3. 
nz  h/r. 

(I  -  Ojjy  ^y*)* 

modulus  of  rigidity  of  the  fscings* 
moduli  of  rigidity  of  the  cor*. 

^xy  <fl  +  t2) 


*xy‘  ffy* 


6  r/h. 


6  r/h. 
o 


Poi.eon' *  ratios  of  the  facing*, 
defined  by  equation  (51), 


S  /  S  . 
y1  °x 
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Figure  3.  -  -Section  of  a  cylindrical  ahell,  where  r  ia  the 
radiua  of  the  middle  aurface  of  the  ahell,  c_  ia  the  thick 
nee  a  of  the  core,  f;  and  f^  are  the  thickneaaea  of  the 

facinga,  £  ia  the  diatancc  indicated  in  the  figure,  and 
£  ia  the  coordinate  indicated  in  the  figure. 

Note:  The  curved  linea  are  area  of  concentric  circles. 


Figure  5.  --N  value*  for  buckling  of  cylinder*  of  sandwich  construction*  in  axial  compression. 
Isotropic  sandwich  facings  of  equal  thickness,  orthotropU  tore. 


e> 

N  !VlN3Wd3dX3 


maammjmsm. 

of  the 

TorMt  product#  Iebaratorj,  forest  Berrioa 
U,  8.  Dapartaant  of  Agriculture 


J.  Aima  h*u . . . 

7«.  J .  HrlwTitt,  c  £  c 

n«rss?e  Boss  Sieffss,  Is  Charts 

Claude  A.  Brown,  Chief.......... 

Gardner  H.  Chideater,  Chief . 

Donald  0.  Colaean,  Chief. ....... 

Herbtrt  0.  Tlaiachar,  Chief . , 

H.  P.  A.  John# on,  Chief . . 

Kenneth  W.  Kruger,  Chief........ 

Ralph  M.  Llndgren,  Chief......... 


. Director 

. . . ....Asfllai&ct  Blfb4tor 

....... .llbnanr  (nadiam  Branch  iSoun ) 

.... ... it... .Fla oal  Control 

. . . . . . . Pulp  and  Paper 

Rea  oaroh  PuMaoationa  and  Itfomtlw 
................... .Tlabar  Prootaalng 

. . . . fhyaloa  and  Engineering 

, . . ........  . Packaging  Reaearoh 

. . . ..Wood  Pregerratlon 


Edvard  0.  hooka,  Chief . . . . . . . Wood  Chanda  try 

Cordon  D,  Lagan,  Chief.. ..................... Adnialatratlre  Haaag— ant 

Harold  L.  Hitohall,  Chief. . .  glebe  r  growth  and  Utilisation  Halation# 


